Theory of Spin-Resolved Auger-Electron Spectroscopy from Ferromagnetic 

3d- Transition Metals 
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CVV Auger electron spectra are calculated for a multi-band Hubbard model including correlations 
among the valence electrons as well as correlations between core and valence electrons. The interest 
is focused on the ferromagnetic 3d-transition metals. The Auger line shape is calculated from a three- 
particle Green function. A realistic one-particle input is taken from tight-binding band-structure 
calculations. Within a diagrammatic approach we can distinguish between the direct correlations 
among those electrons participating in the Auger process and the indirect correlations in the rest 
system. The indirect correlations are treated within second-order perturbation theory for the self- 
energy. The direct correlations are treated using the valence-valence ladder approximation and the 
first-order perturbation theory with respect to valence-valence and core-valence interactions. The 
theory is evaluated numerically for ferromagnetic Ni. We discuss the spin-resolved quasi-particle 
band structure and the Auger spectra and investigate the influence of the core hole. 



I. INTRODUCTION 



Auger-electron spectroscopy (AES) and the comple- 
mentary appearance-potential spectroscopy (APS) have 
become valuable tools for investigatiiM, pthe,|-eiectronic 
structure of solids and solid surfaceairaa'QBQfl'El. They 
represent highly element specific and non-destructive 
methods with a comparatively simple experimental set 
up. The Auger line shape from a core-valence-valence 
(CVV) process yields information on the occupied part 
of the valence band, while the APS provides insight into 
the unoccupied valence states. However, much effort has 
been sjjient on the detailed interpretation of the spectra. 
LanderEI suggested that the spectrum obtained by AES 
(APS) is given as the self-convolution of the occupied 
(unoccupied) valena&, density of states (DOS). On the 
other hand, Powellty discovered the CVV Auger line 
shape of Ag to behave "anomalously" in the sense of 
Lander's self-convolution model. These anomalous fea- 
tures are by now well known to be caused by correla- 
tion effects dominating the electronic properties of vari- 
ous solids. Therefore, AES (APS) seems to be a useful 
technique to study electron-correlation effects, but it is 
doubtful whether it is able to compete with one-particle 
spectroscopies, such as photoemission (inverse photoe- 
mission), in deriving the DOS by deconvolution. 
In the theoretical treatment of the CVV Auger process, 
there are mainly two problems. The first one is to take 
into account the correlation effects. Here one may distin- 
guish between the direct and indirect correlations. The 
direct correlations describe the correlations of those elec- 
trons which participate in the Auger process. They are 
responsible for the most prominent effects in the Auger 
line shape as compared to the self-convolution. On the 
other hand, the indirect correlations among the elec- 
trons in the rest system manifest themselves in the quasi- 
particle density of states (QDOS) as a renormalisation of 
the one-particle DOS. 
The second problem is the calculation of the transition- 



matrix elements for the Auger process as well as the scat- 
tering of the outgoing Auger electron (cf. Refs. 0,12h. 
These effects will (slightly) modify the bare line shape 
and may become important for a refined interpretation 
of experimental data. Within the present paper, however, 
we set aside this second problem and like to concentrate 
on electron-correlation effects in AES from ferromagnetic 
3d-transition metals. 

Within the framework of the single-band Hubbard model, 
correlation effects can be treated exactly for systems with 
completely filled op-empty bands, as was first shown by 
Cini and SawatzkyEjIiJiij. The generalisation to the case 
of degenerate bands was introduced in Ref. Eq and further 
analysed in Ref. |l^, for example. These results may-also 
be extended to include the core-valence interactior£3. 
Considering the more general case of partially filled bands 
introduces several complications concerning the indirect 
as well as the direct correlations. For the indirect valence- 
valence correlations there is a number of approximation 
schemes applicable to a multi-band Hubbard model. A 
method which reproduces the experimentally observed 
Curie temperature quite wdl. especially for Ni, is the 
spectral density approachllSEJ. Other approaches are, 
for example, the generalisation of the single-band npndi-. 
fied perturbation theoryO to the multipband modeEj'E^ 
and Quantum-Monte-Carlo simulatiopaHJ in connection 
with the dynamical mean-field theorycj. However, these 
methods suffer from some necessary restrictions concern- 
ing the completeness of the Coulomttpiatrix. This is not 
the case for the_fluctuation exchangee^ and the Hubbard I 
approximationcj, for example. For a more detailed dis- 
cussion on the indirect valence-valence correlations see 
Ref. ||. 

For the treatment of the direct correlations, one has the 
exact-diagonalisation methndr.1 for small systems and the 
equation-of-motion methodE3c3 with its in general un- 
controlled termination of the hierarchy of the equations 
of motion. Another approximate solution is 
valence ladder (VV ladder) approximatioi 



and its generalisation to include the core-valence inter- 
actionE£(L3EIIca (CVV ladder). In particular one has to 
account for the broken translational symmetry in the ini- 
tial state of AES, caused by the presence of the core 
hole and its screening due to the valence electrons. In 
the final state this interaction is responsible for the sud- 
den response of the valence electrons due to the destruc- 
tion of the core hole. In the limit of completely filled 
or empty bands the ladder approximations recover the 
above-mentioned exact solution. 

Here the interaction strength is taken as the control 
parameter, which is correct in the weak-coupling regime. 
We are aware that this method has restrictions for the 3d- 
transition metals. However, in this work we prefer a com- 
mon treatment of one-particle spectroscopies (photoemis- 
sion and inverse photoemission) and two-particle spectro- 
scopies like AES. To be concrete, we will use the second- 
ordor .perturbation theory around the Hartree-Fock solu- 
tiorOCj'El for the indirect valence- valence correlations. 
The direct correlations will be treated by applying two 
different methods, i. e. the VV ladder approximation and 
the first-order perturbation theory in the valence- valence 
and core-valence interaction. 

Within this approach it is possible to include a real- 
istic one-particle input taken from tight-binding band- 
structure calculations. We do not only account for the 
degeneracy of the 3d-band but also for the hybridisation 
with the 4s and 4p states. The theory is formulated and 
evaluated for a non-orthogonal basis set where the states 
can be distinguished by the angular momentum quan- 
tum number and the cubic harmonic index. This facili- 
tates the interpretation of the resulting spectra. Further- 
more, we do not restrict ourselves to correlations among 
the final-state holes only and include core-hole effects 
from the very beginning. This implies the necessity for a 
proper treatment of the initial state where the core-hole 
screening breaks the translational symmetry. The theory 
is implemented numerically and evaluated for ferromag- 
netic Ni. 

The paper is organised as follows. In the next section 
we will i ntroduce the model under consideration. In sec- 
tion III we give the expression for the Auger intensity. 



Section [V concentrates on the indirect and section ^ on 
the direct correlations. Finally, section VI concludes the 



paper. Some details concerning the non-orthogonal basis 
set are given in the appendix ^ 
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The index i refers to lattice sites, a is the spin-index 
(cr =T,i), L = {/,77i} is the orbital index with angular 
momentum quantum number I and cubic harmonic in- 
dex m. c\^^ {ciLa) denotes the creation (annihilation) 
operator of a valence electron at the lattice site i with 
spin a and orbital index L while b\^ (bia) creates (an- 
nihilates) a core electron. The core states are assumed 
to be non-degenerate and dispersionless with the one- 
particle energy ec well below the chemical potential fi. 
The hopping integrals tf^f 



{iLcj\h^yL'a')=tfJ^'5^, 



(2) 



{h^^ denotes the Hamiltonian of the tight-binding band- 
structure calculation) are taken from Ref. ^ as well as 
the overlap integrals Sf^^f" 



{iL(j\i'L'a') = Stir 5, 



(3) 



if^f and S^^f" refer to a non-orthogonal basis set (see ap- 
pendix H) . Contrary to an orthonormal basis set (where 
the overlap matrix is replaced by Su'Slli), the basis 
states under consideration can be characterised by the 
orbital index L = {l,m}. The construction operators 
likewise refer to the non-orthogonal basis and satisfy the 
following anti-commutation rules: 



[CiLaj^i'L'cr'l^ — U 



CiLa, C^ij^,^ 



= is-')-"' So 



(4) 



It should be noted that the action of the creation operator 
on the vacuum stat e cli^„\^), in general, does not yield 
\iLa) (see Eq. (A. 5) of the appendix^. 

To describe the correlations among the valence elec- 
trons (VV) as well as the correlations between valence 
and core electrons (CV) the interaction consists of two 
parts 
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CV 
dc ■ 



dc 



(5) 



II. MODEL 

The Hamiltonian Ti = i?o ^ l^-N + Hi is decomposed 
into a one-particle part Ho — /iiV and an interaction part 
H\. N is the operator for the particle number. The one- 
particle part describes non-interacting valence and core 
electrons: 



Ho-^iN = J2 



LL' 



nS, 



LL' 



'^iLa'^i'L'iy 



Here the occupation number operator for valence elec- 
trons is riiLa = Cij^^CiLa and for core electrons n'l^ — 
bl^bia-. Assuming a strong screening of the Coulomb in- 
teraction, the interaction part is taken to be purely local. 
UliL^l^Ls are the on-site Coulomb-matrix elements for 
the valence electrons. 

The electronic structure of the 3d-transition metals may 
be understood considering mainly two types of electronic 
orbitals: The 4s- and 4p-states which form broad free- 
electron like bands. They should be well described by 



the band-structure calculation. The other group are the 
well localised 3d-states which in the solid form relatively 
narrow bands positioned around the Fermi energy. The 
localised nature of the 3d-electrons gives rise to impor- 
tant dynamic 3d-3d correlation effects which are believed 
to be responsible e. g. for the magnetic behaviour of the 
3d-transition metals. These correlations may not be ad- 
equately taken into account within a mean-field picture. 
We thus treat them separately. 

Exploiting atomic symmetries, one is able to express all 
remaining Coulomb-matrix elements for th©,3d-electrons 
in terms of three effective Slater integralsEd'cfl (F^, F^, 
F^) only. These integrals are connected to averaged val- 
ues for direct 



U 



sE^- 



25 



LL' 



F" 



L,L' 



and exchange interaction terms 



•^ - ^ E ^^^'^'^ 



L^L' 



F^ + F^ 
14 



(6) 



(7) 



For 3d-elemeni|S-.one has to a good accuracy the ratio 
F^/F^ Ki 0.62Ti to be that of free ionM. U and J are 
treated as free parameters to be fixed by comparison with 
experimental results (see section |lV|) . 
The CV interaction part is necessary to describe the core 
hole effects in AES. t/£ in equation (J^) are the Coulomb- 
matrix elements between the valence and the core elec- 
trons which can be fixed by assuming complete screening 
of the core hole by the valence electrons (see section IV). 
To avoid a double counting of interactions, we subtract 



the correction H 



VV(CV) 
dc 



which is to a good approxima- 



tion jtiie Hartree-Fock part of the respective interaction 
termcj. 



III. AUGER INTENSITY 

The Auger process can be divided into two subpro- 
cesses. The first one is the creation of a core hole with 
spin Uc at the lattice site ic by absorbing an x-ray quan- 
tum for example. The second subprocess is the radiation- 
less decay of the core hole via ejecting an Auger electron 
with spin a and momentum k. Provided that the life 
time of the core hole is large compared to typical relax- 
ation times of the valence electrons in the presence of the 
core hole, the two subprocesses become independent from 
each othem (two-step model). This implies the absence of 
any decay term in the Hamiltonian. Consequently n1 ^ 
is a good quantum number, \H, '^1 a\^ =0. We can 
concentrate on the second subprocess. Within the two- 
step model the initial state for the Auger transition pro- 
cess is the ground state within the subspace 'W of the 
Hilbert space Ti that is built up by all many-body states 
with n1 ^ = 0. To perform thermodynamic averages in 
practice, one has to take into account this restriction by 
introducing an additional Lagrange parameter. 



The transition process itself is described by the tran- 
sition operatoiO 



^ k(T(T, 



.-E^i< 



-f-2 t t !^. 



(8) 



L\,L2 



k and a denote the quantum-numbers of the Auger elec- 
tron, (Tc the spin of the core state involved. The (intra- 
atomic) Auger-matrix elements are given by 



M/ 



(X 



(icLi,ici2|-ffcoulomb|k, ic) (9) 

1 



|ri - V2 



■*k(ri)0(r2-R,;J 



where ^ is the valence orbital, <& the one-particle wave- 
function of the Auger-electron and (f) the core state. The 
bar denotes complex conjugation. 

Following Ref. BS we consider the (retarded) three- 
particle Green function, relevant for AES, which is de- 
fined as (with the abbreviation L = (Li,L2) and L' = 



G'iLXE) = {{Tl.;T^ 



ko"0"c 



J)i 



'kcra,,' -'kcrcr. 



(10) 



E^tkG5'"">)^^^k 



L,L' 



with 



((.; .))e refers to Zubarev Green functiongla'ta. The 
AES intensity is then mainly given by the three-particle 



l(3) 



(lA)3G|fi.(£^): 



spectral-density Aj^j^ (E) = 

haa^ {E + e,-fi)<^ 5{E - E{'k))A'il^^ (E) . (12) 

Here -E(k) is the dispersion of the Auger-electron. 

In general the three-particle Green function will be a 
(complicated) functional of one-particle Green functions. 
This functional represents the direct correlations. In 
the following we concentrate on the indirect correlations 
first, i. e. on the determination of the relevant one-particle 
Green functions. 

IV. INDIRECT CORRELATIONS 
A. Valence-band interaction 

We consider the (retarded) one-particle valence-band 

Green-function {{ciLa] Ciii^,„))E- Using a matrix notation 
with respect to the orbital index L = {I, m} 



(x,,o 



LL' 



yLL' 



(13) 



and defining a lattice-Fourier transformation 



Xk 



wX' 



-?k(Ri-R.^, 



^ii' 1 



(14) 



where iVg is the number of lattice sites, we get Dyson's 
equation in the form 



(ji. 



1 



■W = 7^E 



„ik(R,-R;, 



(S + M)Sk-tk-Sk^(-B) 



(15) 



Here the matrices tk and Sk are the Fourier-transformed 
hopping and overlap integrals of Eqs. (|^) and (||), re- 
spectively. SkCT(-E) is the self-energy. The one-particle 
spectral density is given by 



^■n'a{E) 



1 



^Gu'a{E + iO-^ 



(16) 



The on-site terms of the spectral density are diagonal in 
the orbital index as a consequence of lattice symmetries, 
and we have for the orbital resolved QDOS 



p'^AE) ^ A^^AE ~ ^i) 



(17) 



where we dropped the site index. The total QDOS is 
obtained via 



Pa{E) 



/^ 



Ti{S,,A,,,,(E~ii)} . 



(18) 



To calculate the self-energy, we use a standard ap- 
proximation scheme, the second ordp,|P4ptuthajtion the- 
ory around the Hartree-Fock solutionEEaJid'ElEl (SOPT- 
HF). It is known from the single-band Hubbard model 
that the non-local terms of the SOPT-HF self-energy 
rapidly decrease, with increasing number of shells taken 
into accountOO. Furthermore, due to the band degen- 
eracy there is a much weaker k-dependence of the jSOPT- 
HF self-energy compared to the single-band casel3. We 
may thus employ the local approximation 



^LL' 



Eilf^JE) « I],^> {E) Sw = KiE) 5Luh. 



(19) 



As for the on-site Green function, lattice symmetries re- 
quire the on-site self-energy to be diagonal in the orbital 
index. 

The Hartree-Fock (HF) contribution to the self-energy 
reads as 



yj(HF),L 



(0)Li 



Y.{uLL.LL\n^J^-n^^r)+ (20) 

{Ull.ll, - Ull,l,l) (n^^ - 4"^^^) } 



where n^ = {niLa) denotes the expectation value in 
the full model and nJr the expectation value of the 
band-structure calculation that stems from the double 
counting correction in equation (p[). Approximating the 
self-energy by E^ , corresponds to the LDA-|-[/ ap- 
proached. Here we additionally include the next order in 



the interaction. The second-order contribution (SOC) to 
the self-energy reads as 

dx dy dz 



-^(SOC),L 



(E)^ 



(21) 



E — X + y — z 
X (/_(a;)/_(-jy)/_(z) + /_(-a;)/_(y)/_(-z)) 

Y^ Ull,LsL2 Pa^{x) [Ul,l2LL, p^^„{y)p^„{z) 



L\.L2;Li 



with the Fermi function f-{E) = {e^^ + 1)^^ and where 

the local HF spectral-density p^{E) = A^f^'^^{E) is ob- 
tained by using the HF self-energy ( pO| ) in equation ^a). 
The SOPT-HF self-energy 

S,(£;) = SSH^) + S(foc)(^) (22) 

determines the full Green function via equation (|n 



B. Core- valence interaction 

1. Valence-band Green function 

Let us now focus on the core hole screening in the ini- 
tial state for AES. The CV interaction and the presence 
of the core hole introduce an additional (Hartree-like) 
term 



,(CV,c),L _ 



= -SuMf 



(23) 



to the valence-band self-energy which breaks the trans- 
lational symmetry. This term represents the core-hole 
potential at the lattice site ic where the core hole was 
created (the superscript "e" indicates averaging in the 
restricted Hilbert-space H"^; see section IH). It is respon- 
sible for the screening. The valence-band self-energy in 
the presence of the core hole then reads 



^UE) 



^(VV,o) 



(E) 



,(CV,e) 



(24) 



,(VV.o). 



Sct ' {E) incorporates the VV-correlation effects and 
has the same structure as the self-energy (^ for the 
translational invariant system. But in contrast to its 
translational invariant counterpart, SJt '" (E) is defined 
in terms of the Green functions in the presence of the core 
hole. 

The valence band Green function G°j,g.(i?) in the pres- 
ence of the core hole can be obtained by using Dyson's 
equation in the form 

Gl,AE) = Gu'AE) (25) 

+ Y,G,ME) {^AE) - ^a{E)) G],AE) . 
j 

In general 'E^^{E) — 11cr{E) ^ for a certain num- 
ber of shells around the core-hole site because the VV- 
correlation effects depend on the occupation numbers, 



which as a consequence of the screening locally differ from 
the translational invariant ones. In the following we as- 
sume complete screening, i. e. charge neutrality at the site 
«c, which is reasonable especially for 3d-transition metals 
because the screening time scale is small compared to the 
life-time of the core holeEj. This impHes S^^ (E) - T,a (E) 
to be small for all sites i except for i = ic- Neglecting 
the terms for i ^ ic one can solve equation (|2£ 



G\,„{E)^G,,,{E)^G. 



it^a{E) 
-1 



[i:i,{E)-i:,{E)y'-G,^,ME) 



G, 



(26) 
'4E). 



For i = i' = ic one obtains the local screened Green- 
function at the core-hole site: 



G^ 



r{E) = 



1 



G-XJE)-{^1^{E)~^,{E)) 



(27) 



The assumption of complete screening will be utilised as 
condition to fix the CV-interaction parameter, which is 
taken to be the same for s-, p- and d-orbitals (C/£ = Uc). 



> 







-4-2 2 -6 -4 -2 

(E -ii) (eV) (E -ji) (eV) 



FIG. 1: Spin-resolved QDOS per atom for s-, p-, t2g-, eg- 
states and total QDOS for U = 2.47 eV, J = 0.5 eV and 
T = OK. Left: unscreened QDOS.rThin dotted line: tight- 
binding band-structure calculations for paramagnetic Ni. 
Right: screened QDOS at the site ic in the presence of the 
core hole, Uc ~ 1.81 eV. 



2. Core Green function 

To take into account the CV interaction for the core 
Green function 



1 



(28) 



E + ^I-ec-J:l,^iE) 



one may calculate the core self-energy S^^^^ {E) using e. g. 
the SOPT-HF in the same way as for the valence-band 
Green function. On the other hand it is believed that the 
core states are influenced by other and presumably more 
important effects, such as life-time effectsB. In fact, the 
core spectral density obtained within SOPT-HF turns 
out to be dominated by a delta-peak that is shifted by 
about 1 eV below Ec — /i. This does not affect the Auger 
line shape. Therefore, we assume for convenience the core 
self-energy to be zero. The spectral density becomes 



a»c<Te(^) 



1 



<ig,^,^{E + iO+) ^ S{E + n ~ ec) . (29) 



C. Results for Ni 



Before we discuss the results for fcc-Ni, we like to make 
a short remark concerning the numerical evaluation of 
the theory. The k-sum for the local Green function in 
equation ( |l5| ) was performed on a mesh of 240 k-points 
within the irreducible piart of the Brillouin zone using 
the tetrahedron methocO generalised to complex band 
structures, similar to that presented in Ref. IsO. The eval- 
uation of the total QDOS (|l|) as well as the QDOS in 
the presence of the core hole was done in k-space. For 



the latter, equation ( pSf) has to be used to perform the 
Fourier transformation. 

The effective Slater integrals or, equivalently, the aver- 
aged direct and exchange interaction parameters are cho- 
sen a.s U — 2.47 eV and J = 0.5 eV. This leads to a cal- 
culated magnetic moment per atom of tti = 0.56/iB|-ajb 
T — OK which is the same as the measured momentEiJ. 
With the ratio J/[/ « 0.2 we assume a typical value for 
the late 3d transition metals. The values givenjin the 
literature, for instance II, = 3.7 eV, J = 0.27eVli2l and 
U = 2.97 eV, J = 0.8eVc3, are of the same order of mag- 
nitude but slightly overestimate the magnetic moment 
within the present theory. 

The "free" DOS, used as starting point for our theory is 
shown on the left of figure H (thin dotted line) and corre- 
sponds to tight binding band structure calculationscj for 
paramagnetic fcc-Ni. 

The left-hand side of figure |l] shows the QDOS per 
atom for the model parameters given above. As is known 
from the experiment, Ni is a strong ferromagnet, i. e. the 
majority-spin states are fully occupied. The renormali- 
sation effects of the a priori uncorrelated s- and p-states 
seen in figure ll| can be traced back to the hybridisation 
with the d-states. 

Taking into account the presence of the core hole and 
following the procedure to fix the CV interaction param- 
eter as described above (charge neutrality), leads to Uc — 
1.81 eV. The corresponding "screened" QDOS is plotted 
on the right-hand side of figure |l|. The structure has 
changed remarkably. Spectral weight of the d-electrons 
from the upper band edge is transferred to lower ener- 
gies, the s- and p-states are more populated, too. The 
screened magnetic moment at the site ic {ml = 0.1/ib 
at T = K) is considerably decreased since the local oc- 



cupation is increased (see also figure ^. 

The left side of figure g shows the self-energy E^ (E) . 
Within an energy range of about 1 eV above and be- 
low E ^ OeV one has 3S^(£') oc E^. Thus we have 
well-defined quasi-particles at the Fermi energy and their 
weight 



dm:^{E = 0) 



dE 



(30) 



is 0.887 for the t2g-t-states and 0.893 for t2g-i-states. For 
eg-states we find 0.878 and 0.883, respectively. For ener- 
gies above — 2eV, where one finds clearly distinguishable 
structures, a significant band-narrowing caused by the 
real part of the self-energy is observed. While the imag- 
inary part of the self-energy leads to strong damping ef- 
fects in the QDOS (figure |l|, 1. h. s.) for energies below 
— 2eV. About — 6eV below the Feoai energy where one 
expects the "Ni-6-eV-satellite"E3aO we find the largest 
damping effects. However, we do not find the correlation- 
induced 6 eV satellite. This is not surprising by apply- 
ing the SOPT-IIELpi any other finite-order perturba- 
tional approachEHlLSEZl. For different interaction parame- 
ters (larger [/, smaller J) a small shoulder in the QDOS 
of the d-states is visible as was also mentioned in Ref. ^. 
On the right-hand side of figure the self-energy in the 
presence of the core hole Yf^^^[E) is plotted. Again there 
are well-defined quasi-particles, but with an enhanced 
weight compared to the case where the core hole is absent 
(<^f = 0.940, zZl = 0.934, ^^^ = 0.937, z°;^ = 0.936). 
The screened case behaves less correlated than the un- 
screened case since here one is closer to the limit of com- 
pletely filled bands. 

Finally, we show the local magnetic moment per atom 
as a function of temperature in figure 0. The mag- 
netisation curves (figure Q, 1. h. s.) have a Brillouin- 
function-like form, except for the eg-magnetisation which 
shows up a maximum at T « HOOK. This can be 
traced back to a transfer of charge carriers from the eg- 
orbitals to the t2g-orbitals with increasing temperature. 
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Because Ni is a strong ferromagnet the charge-carrier 
transfer leads to an increase of the Cg-magnetisation. 
Contrary, the t2g-magnetisation is decreased in addition 
to the usual temperature-induced depolarisation. This 
leads to a temperature-dependent increase of the ratio 
"^og/(wt2g +"T-Cg) as itisJmown from polarised neutron- 
scattering experimentsE§E3. For T = OK this ratio is 
0.20p 
0.19E: 



in good agreement with the measured value of 
As is observed experimentally, the s- aiid.ai- 
states couple antiferromagnetically to the d-statesESEj. 
The Curie temperature turns out to be Tc = 1655 K and 
is thereby about a factor 2.6 larger than the measured 
value of 624 KE3. The large value for Tc is probably due 
to the mean-field character of the SOPT-HF. Note, how- 
ever, that a simple LDA+C/ (HF) calculation for the same 
parameters U and J yields a T = magnetisation 0.57 /.te 
and a Curie temperature of approximately 2500 K. 
The local moment at the site ic as function of temper- 
ature is shown on the right-hand side of figure |3[ Its 
d-contribution is strongly reduced compared to the un- 
screened case while the s- and p-moment is increased. 
The strong reduction of the total magnetic moment 
will influence the spin polarisation of AES, since the 
"screened" QDOS enters the Auger Green function [lC| . 
Note that the total magnetisation has to be calculated 
using equation (HH) , incorporating hybridisation with de- 
localised states. 



V. DIRECT CORRELATIONS 

To express the Auger intensity as a functional of 
the one-particle Green functions we consider the dia- 
grammatic expansion of the three-particle Green func- 
tion Gf^^l^'{E) (L = {Li,L2)). Here we can re- 
strict ourselves to the direct diagrams and incorporate 
the exchange diagrams by introducing "direct" {D\'-^^ = 

and "exchange" Auger matrix elements {E^-^^ — 
The Auger intensity then reads as 






Ikaa, {E + €c- tJ.) (^ 



SiE-Eik)) 



L,L' 



D, 



(31) 
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FIG. 2: Real and imaginary part of the self-energy. Left: for 
the translational invariant system (Ect(_B)). Right: for the 
system in the presence of the core hole {'S1^^{E)) 



FIG. 3: Left: magnetisation as function of temperature. 
Right: local moment as function of temperature in the pres- 
ence of the core hole. 



The Auger matrix elements are taken to be constant. 
Following Ref. ^we set: 



M 



L1L2 



1 for Li < L2 
— 1 for Li > L2 



(32) 



We thereby account for the singlet contributions, i. e. the 
holes in the final state have opposite spin {a = — Cc), as 
well as for the triplet contributions {a — ac) to the Auger 
intensity. The triplet contributions would be ignored if 
the Auger matrix elements were chosen to be symmet- 
ric in the orbital index (M^j^^= = M^^^^) because the 
transition operator then vanishes, as can be seen from 
equation (||). 



A. VV ladder approximation 



convolution of the partial QDOS 

At'^^'-'iE) ^ 5^^, jdE'pl\{E-E')pllSE') (35) 
X (/_(£' - E)f^{-E') -UE- E')UE')) 
using the spectral representation: 

In equation (p4) we have applied the local approxima- 
tion, i. e. only the on-site elements for i ~ ic are assumed 
to be non-zero. This approximation is analogous to the 
local approximation for the one-particle Green function 
in section 0. 



We consider two different approximation schemes for 
the treatment of the direct correlations. In the first 
approach, following Refs. pO|j3l|j3^ , P3| , p4] we neglect the 
direct CV correlations and treat the direct VV cor- 
relations by means of the VV ladder approximation, 
which becomes exact in the limit of completely filled or 
empty bands. A typical diagram contributing to the 
VV ladder is shown in figure 0. The solid lines rep- 
resent the renormalised one-particle propagators of the 
valence-band while the wiggly line is the one-particle 
core-propagator. The dashed line corresponds to the VV 
interaction. Summing up all diagrams of this kind, yields 
the VV ladder approximation. The three-particle spec- 
tral density is given by 

^d^{E) = JdE'A^X^'iE + E') 

xa^^aAE')f+iE + E') (33) 

=^ A^J^^^'{E + e,-p)UiE + e,-^^) 

where f+{E) = (e^^ — 1)^^ is the Bose function. The 

two-particle valence-band spectral density A] 'J^ (E) 
is obtained from the corresponding two-particle Green 
function. Using a matrix notation with respect to L = 
{Li,L2} and L', the two-particle Green function reads 
as 

g5.. (E) = g5"J^ (E) (1 - Ug5"J^ (E)) "' (34) 

with C/ll' = UliL2L{L'^- Cfj^'^^ {E) is the two-particle 
Green function which has to be calculated form the self- 



B. VV and CV correlations 

A straightforward way to include the direct CV cor- 
relation on the same level as the direct VV correlations 
has been discussed in Refs. |35| , p6| , p7|j3q . This leads to 
the ("three-particle") CVV ladder approximation. For 
the limiting case of completely filled or empty bands the 
CVV ladder represents the exact solution and recovers 
the VV ladder but— shifted energetically by 2Uc due to 
the CV interactionta. A typical diagram is shown in fig- 
ured where the dotted line represents the CV interaction. 
The CVV ladder approximation leads to a coupled set of 
Fredholm integral equational3. For a multi-band model, 
however, the numerical evaluation is beyond our present 
computational capacities. We therefore discuss a simpler 
approximation where only diagrams up to first order in 
the direct correlations are retained (see figure 1^) . 



C. Results for Ni 

The calculated Auger spectra for Ni resulting from dif- 
ferent approximations are shown in figure uv The. core 
hole is assumed to be unpolarised (non-resonantEil pro- 
cess). The intensities for core spin Uc and — CTc are added 
incoherently: h{E) = {I^^^{E) + I,_,^{E))/2. How- 
ever, the Auger intensity is still spin-dependent due to 
the ferromagnetic order in Ni. 

In figure ^ we plotted the total Auger-intensity /| {E) + 
Il{E) on the left and the spin asymmetry {Ii^{E) — 
Il{E))/{I^{E) + IiiE)) on the right. Part (a) shows the 
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FIG. 4: Typical diagram of the VV ladder. Solid line: 
renormalised valence-band propagator. Wiggly line: core- 
propagator. Dashed line: VV interaction. 



FIG. 5: Typical diagram in the CVV ladder approximation. 
The notation is the same as in figure M. The additional CV 
interaction is represented by the dotted line. 



FIG. 6: Diagrams up to first order in the VV and CV inter- 
action (notation as in figure q). 



result of the self-convolution model (equations (33) and 
( P5| ) inserted in (|3l|)), i.e. the self-convolution of the oc- 
cupied QDOS in figure 13 (1. h. s.). Direct correlations 
and core-hole screening are neglected altogether. Part 
(b) corresponds to the VV ladder approximation starting 
from the (unscreened) QDOS. Taking additionally into 
account the screening effects introduced by the presence 
of the core hole in the initial state results in (c). The 
spectrum obtained by the first order in the direct VV 
and CV correlations (see figure o) and by the screened 
QDOS is plotted in part (d). 

As one can see in the plots on the left-hand side, the VV 
interaction is too weak to produce bound states, no sharp 
satellite appears, and the spectra appear to be band-like. 
Compared to the self-convolution (a), however, a consid- 
erable shift to lower energies is observed in (b) . This shift 
results from the direct correlations between the two final- 
state holes in the valence band. In (c) the main peak is 
shifted to still lower energies. This is an effect of the core- 
hole screening in the initial state and can be traced back 
to the redistribution of spectral weight in the screened 
QDOS (figure |l]). Compared to (a) and (b), the total 
AES intensity is clearly increased in (c). Again this is a 
consequence of the core-hole screening since the number 
of occupied states available for the Auger process is in- 
creased (figure n^, r. h. s.). The spectrum shown in (d) not 
only includes the initial-state core-hole screening but the 
final-state effects due to the destruction of the core hole. 
Compared to (c), where the initial state is described in 
the same way, these effects result in a strong shift of the 
main peak to higher energies. This shift almost exactly 
compensates the shifts to lower energies that are due to 
direct VV correlations (b) and the core-hole screening 
(c). However, a weak shoulder at about E = — 8.5eV 
remains in the spectrum (d). 

In all cases there is a high spin asymmetry (up to —50 %) 
for energies between approximately — O.SeV and 0. This 
is a consequence of the fact that Ni is a strong ferromag- 
net. There are almost no t-electrons above about — 0.4 eV 
(see figure |l|) that can participate in the Auger-process. 
The main contribution to the intensity is therefore due to 
triplet configurations where the two final-state holes or, 
equivalently, core- hole and Auger electron have spin J,. 
However, the intensity is very small in this energy re- 
gion. By taking into account the screening of the core 
hole in the initial state (compare (b) and (c)) the spin 
asymmetry is reduced over the whole energy range which 
essentially is the same effect as the reduction of the lo- 
cal magnetic moment at ic caused by the presence of the 
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FIG. 7: Left: total intensities It{E) + Iy{E). Right: 
spin-asymmetry (/t(-E) - Ii{E))/{I^{E) + Ii{E)). (a) self- 
convolution without screening of the core hole in the initial 
state; (b) VV ladder without core hole screening; (c) screened 
VV ladder; (d) direct VV and CV correlations included up to 
first order (figure 0). 



core hole (figure 0) . The total spin-polarisation 



P 



JdEmE)-JdEI^iE) 
JdEmE)+JdEI^{E) 



(37) 



in the case (d) is 2.6% and 1.6% foe— (c). Both val- 
ues are close to the experimental valua£^ of 2 % for the 
MiM45M45-process. The cases (a) and (b) with a polar- 
isation of 8.7% and 9.3%, respectively, overestimate the 
total spin-polarisation compared with the experimental 
value. 

For the calculation of the orbitally resolved contribu- 
tions to the Auger intensity we may restrict the summa- 
tion in equation ( p3l| ) to orbital indices (Li, L2, L'l, L2) 
belonging to t2g (CgJ^character only. The resulting con- 
tributions are shown on the left (right) of figure 0. The 
contributions due to the remaining terms are plotted in 
the middle. 

In all cases (a)-(b) the t2g contributions are clearly 
stronger compared with the Cg contributions. The ratio 
between the t2g and Cg partial intensities corresponds to 
the different degeneracies. Comparing the cases (a)-(d), 
we notice that there are essentially the same trends in 
the partial intensities as for the total intensities, and the 
discussion is the same as above. The line shape in case 
(d), however, shows up some fine structure, especially in 
the t2g partial intensity, which is not that pronounced in 
the total intensity. The shoulder at the low-energy tail 
of (d) is due to direct VV correlations and may be inter- 
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FIG. 8: Contributions to the total Auger intensity (figure M) of 
processes involving t2g-electrons only (left), eg-electrons only 
(right) and both kinds of d-electrons (middle). 



preted as a hint for the formation of a bound state of the 
final holes. Except for this shoulder, a surprising simi- 
larity between (d) and (a) is noticed, even for the orbital 
resolved spectra. This might also be due to the small 
number of diagrams taken into account. However, the 
cancellation of effects according to different interactions 
(VV and CV) was also pointed out in Ref. E^. 



VI. SUMMARY 

In this paper we have investigated electron-correlation 
effects on the Auger line shape of Ni as an example 
of a ferromagnetic 3d-transition metal. The starting 
point is a realistic set of hopping and overlap param- 
eters taken from tight-binding band-structure calcula- 
tions. We additionally consider a strongly screened on- 
site Coulomb interaction between the rather localised 3d 
electrons. The respective Coulomb-matrix elements are 
expressed in terms of effective Slater integrals. Choos- 
ing a non-orthogonal basis, a distinction between the dif- 
ferent angular-momentum characters of the valence or- 
bitals is possible. This is necessary for the precise defini- 
tion of the Coulomb-interaction part in the (multi-band 
Hubbard) Hamiltonian and also facilitates the interpre- 
tation of the Auger spectra. Furthermore, we account for 
the core- valence interaction which is responsible for the 
screening of the core hole in the initial state of AES and 
for the sudden response of the valence electrons due to 
the destruction of the core hole in the final state. 

Within a diagrammatic approach, the indirect and the 
direct correlations can be studied separately. The in- 
direct correlations have been treated by second-order 
perturbation theory around the Hartree-Fock solution 



(SOPT-HF). The VV interaction parameters are fixed 
by assuming a ratio J/U ~ 0.2 and by fitting the ex- 
perimentally observed magnetic moment for T = OK, 
leading to U = 2.47 eV and J — 0.5 eV. This is equiv- 
alent to an intra-orbital interaction of Ullll = 3.07 eV 
(L = {2,7712}). The resulting Curie temperature within 
the theory presented here is by a factor 2.6 larger than 
the experimentally observed one but considerably lower 
than the LDA+U (Hartree-Fock) value. 
The core- valence interaction leads to a breakdown of the 
translational invariance in the initial state for AES. The 
interaction parameter Uc — 1.81 eV is fixed be requiring 
charge neutrality at the site ic where the core hole is cre- 
ated. The screening of the additional core-hole potential 
causes a transfer of spectral weight below the Fermi en- 
ergy and thus a considerable reduction of the local mag- 
netic moment. However, the local magnetic moment is 
finite since the s- and p-electrons also contribute to the 
screening. 

To study the direct correlations we have used two dif- 
ferent approaches. The first is the VV ladder approx- 
imation, which results in a band-like Auger spectrum 
with a single maximum. In a second approach we have 
summed up the first-order diagrams with respect to the 
VV and CV interaction. The resulting spectrum shows 
up a shoulder at the lower tail due to the VV interaction. 
Otherwise the line shape is very similar to that obtained 
by the self-convolution of the unscreened QDOS. 
As far as concerns the line shape, we conclude that the 
different correlation effects, VV and CV correlations in 
the initial and the final state, nearly cancel. However, 
a strong effect of electron correlations has been found in 
the orbitally resolved partial intensities and particularly 
for the spin-asymmetry. 

The calculated spin polarisation is in a good ag|Cfiement 
with the measured one of the AfiM45M45-processLj. This 
corresponds to the excitation of a not too deep lying core 
level, i. e. the two-step model should be applicable. The 
reduced but finite spin polarisation of the M1M45M45- 
process (as compared to the band polarisation) can be 
explained by effects of core-hole screening rather than by 
a core-hole polarisation, caused by a resaaant excitation 
of the core electron into the valence bandtil. Future work 
will show whether these findings also apply to other 3d- 
transition metals. 
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APPENDIX A: NON-ORTHOGONAL BASIS-SET 

There are several advantages for using the non- 
orthogonal basis set {\iLcr)}. The Slater-Koster pa- 
ranietersE3 for the two-center approximation (used here) 
are much more accurate for a naa.-orthogonal basis set 
compared to an orthonormal onec3. Secondly, the non- 
orthogonal basis (LCAO basis) is built up from quasi- 
atomic orbitals. One therefore knows the behaviour of 
the basis states under symmetry operations belonging to 
the 0/i-group which eventually results in the fact that 
local quantities, e. g. the on-site Green function, are di- 
agonal in the orbital index. Furthermore, the Coulomb 
matrix elements can be calculated in a highly symmetric 
way by using 3j-symbols in combination with a transfor- 
mation from spheric to cubic harmonics. The unknown 
radial parts of the basis are parametrised by the effective 
Slater integrals (F", F^, F^). Finally, one can ensure 
that the Coulomb interaction acts between 3d-electrons 
only. 

For the formalism of second quantisation, for the 
many-body and Green-function theory, for the proof of 
Wick's theorem and thus the development of the dia- 
gram technique, however, the use of an orthonormal set 
of one-particle basis states is necessary. Therefore, it is 
convenient to derive all expressions that refer to the non- 
orthogonal basis. 



Sal3 : 



(A.l) 



by a (non-unitary) Lowdin transformationcJ from a re- 
lated set of orthogonal one-particle basis states: 



Dcf. 



El 





)S 



-1/2 
f3a ■ 



(A.2) 



The overlap matrix S is Hermitian, {|5)} indeed repre- 
sents an orthonormal and complete basis set. The com- 
pleteness relation can be written as 

ct a/37 

= El«)'^oi(/5|. (A.3) 

Annihilation (and creation) operators referring to the 
non-orthogonal basis may be defined as: 



Dcf. 



^,-U. 



p c^- 



(A.4) 



It is instructive to see how the creation operator acts on 
the vacuum state: 



410) = Y.^WJ"\^') = Y.\P)sTJ" 

/3 /3 

Furthermore, one gets from the transformations ( |A.2| ) 
and ([A!^): 



Y^\a)Z^^Y. \P)S-p'J^S%^c,^Y.\^)c^. (A.6) 

a a/37 ^ 



Thus, an operator in second quantisation has the 
same structure for both, the orthonormal and the non- 
orthogonal basis set. A one-particle operator O, for ex- 
ample, reads as: 



O = y^4 {a\o\a')ca' 

a. a' 
(A.6) ^-^ t / I I '\ 



(A.7) 



The non-orthogonal Green functions are defined as in 
the orthonormal case, e.g. Gaa'{E) = ((cq;^/)) for the 
one-particle Green function. For example, using the non- 
orthogonal version of the fundamental anti-commutation 
rules ([cq,c]j,]+ = S~l^,, see also Eq. (||)) and the equa- 
tion of motion, the non-interacting Green function turns 
out to be: 



G\E)^{{E + ^,)S-T)-^ 



(A., 



where we used the matrix notation {G)aa' — Gaa' etc. 
For the interacting Green function one has (compare with 
Eq. dll)): 



G(F) = ((F + m)S-T-S(F))" 



(A.9) 



In the same way as for the examples given, one may 
use Wick's theorem, develop the diagram technique etc. 
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